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Abstract 



The effect of gravitational radiation reaction on circular orbits around a 

spinning (Kerr) black hole is computed to leading order in S (the magnitude 

of the spin angular momentum of the hole) and in the strength of gravity 

M/r (where M is the mass of the black hole, r is the orbital radius, and 

G = c = 1). The radiation reaction makes the orbit shrink but leaves it 

circular, and drives the orbital plane very slowly toward antialignment with 

the spin of the hole: tan(t/2) = tan(t /2)[l + (61/72) (S/M 2 ) (M/r) 3 / 2 ], where 

t is the angle between the normal to the orbital plane and the spin direction, 

and to is the initial value of t, when r is very large. 
PACS numbers: 04.25.Nx, 04.30.Db 
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The earth-based LIGO/VIRGO network of gravitational wave detectors (which is now 
under construction) will be used to search for and study the gravitational waves from "par- 
ticles" , such as neutron stars and small black holes, spiraling into massive black holes (mass 
M up to ~ 3OOM ); and ESA's planned space-based LISA [|I|] interferometer will do the 
same for inspiral into supermassive black holes (M up to ~ 10 7 M Q ). To search for the inspi- 
ral waves and extract the information they carry will require templates based on theoretical 
calculations of the emitted waveforms; and to compute the waveforms requires a detailed 
understanding of how radiation reaction influences the orbital evolution. 

For several years a stumbling block has impeded computations of the evolution, when the 
orbital plane of the particle is inclined to the equatorial plane of a spinning hole: No practical 
method has been developed to deduce how radiation reaction influences the evolution of 
the orbit's "Carter constant" [@,[|, which governs the orbital shape and inclination angle. 
This Letter describes the first progress on this problem: a "post-Newtonian" gravitational 
radiation reaction force is used to compute the full orbital evolution to first order in S, the 
magnitude of the spin angular momentum of the black hole, and leading order in the strength 
of gravity M/r at the orbital radius r. (Here and throughout, units with G = c = 1 are used.) 
The analysis is restricted to orbits that initially are "circular" (more precisely, orbits which 
have constant radius r — these orbits are circular in the "orbital plane" discussed below, but 
this plane precesses). However, the method can readily be extended to noncircular orbits 0] 
and (with considerably more difficulty) should be extendible to the fully relativistic regime 
r ~ M. 

The computation of the evolution presented here proceeds as follows: First, in the absence 
of radiation reaction, the orbital motion and the associated constants of motion are reviewed. 
Then, the leading order radiation reaction accelerations that act on the orbiting particle and 
on the hole are derived and used to compute the radiation-reaction-induced evolution of the 
constants of motion. Finally, the evolution of the orbit — its shape and inclination angle — is 
obtained. 

The leading order effect of the spin on the (otherwise Newtonian) orbit was deduced long 
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ago by Lense and Thirring j|] (reviewed by Landau and Lifshitz [§]) — though, of course, they 
regarded the central body as a star rather than a black hole. In fact, our analysis does not 
require the body to be a black hole, but since this is the primary case of physical interest, 
the discussion is phrased in terms of a black hole. 

Let spherical polar coordinates, r, 9, and 0, centered on the black hole, be used to 
describe the location of the particle (these coordinates describe the relative separation of 
the two bodies), with the hole's spin along the polar axis. The Lagrangian for the motion 
of the particle (which, for now, does not have to be circular) is given, to linear order in S 
but otherwise in solely Newtonian theory, by 

C = H (r 2 + r 2 9 2 + r 2 sin 2 (£)0 2 ) + ^ - (1) 

where \i is the mass of the particle. In general, an overdot represents d/dt. The entire 
analysis is to leading order in /i. To leading order in S and in M/r, the motion resulting 
from this Lagrangian is the same as in the Kerr metric, which describes the gravitational 
field of a spinning hole. The use of flat space coordinates, which ignores M/r corrections, is 
adequate to leading order. 

Following standard procedure ||, the Hamilton- Jacobi equation associated with the La- 
grangian (P can be shown to have a separation-of- variables (t, <fi, r, 9) solution, which reveals 
three constants of motion: 

E = H (r 2 + r 2 9 2 + r 2 sin 2 (#)0 2 ) - (2a) 

2 T' 

9 „ . „. • 2uS sin 2 9 , , . 

L z = fir 2 sin 2 (6U - — , 2b 

r 

Q + L\ = /A 4 (0 2 + sm 2 (9)<j) 2 ) - Afi 2 Sr sin 2 (#)0. (2c) 

Comparing these expressions to Eqs. (33.31) of Ref. 0, it is clear that E, L z , and Q 
correspond, to leading order in S and in M/r, to the constants of motion for a test particle 
in the Kerr metric: the energy minus the test particle mass, the ^-component (component 
along the spin axis of the hole) of angular momentum, and the Carter constant, respectively. 
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The analysis of radiation reaction, below, is in Cartesian coordinates, x± = r sin 9 cos 0, 
Xi = r sin 9 sin 0, and x% = r cos 9. In these coordinates, the constants of motion (0) become 
(repeated indices are summed) 

£ = 1^-7^. Pa) 

. 2fiS {x\ + sf) 
^ = fie S j k XjXk -3/2 — , (<3 b ) 

r\ _|_ t2 _ 2 • • _ SefykXjXk . . 

t £ijk% j-^k^ilm-^l^m " \ 1/2 ' V / 

(X m X m ) 

Let the orbital plane have inclination angle t (restricted to < 6 < 7r, where t > 7r/2 
corresponds to an orbit counter-rotating relative to the spin), defined as 



cos £ 



(4) 



(Q + ^) 1/2 ' 

The constants of motion admit orbits of constant radius, just as for the Kerr metric in 
Boyer-Lindquist coordinates. One of the Euler-Lagrange equations implies that, for r to be 
constant, dC/dr = 0, which leads to 



r = Mv~ 



1 + 6(S/M 2 )v 3 cosl , (5) 



'Q + LIY 1/2 



(6) 



V /i 2 M 2 J 

The positive root in any square-root is always chosen. 

The following relationship between the constants of motion for circular orbits is easily 
derived: 

E = -i/it» 2 [l - 4(S/M 2 )v 3 cos l] . (7) 

The other two Euler-Lagrange equations predict that 9(f) and 4>(t) are the same as for 
a circular orbit in the case when S = 0, except that <f) is altered to <fi = <p\s=o + 25*/r 3 . By 
transforming these angular motions to Cartesian coordinates, the following equations for the 
orbit are obtained: 



x\ — r 



%2 = r 



cos 



, n \ {2St\ . frx , . /2St\ 

(ilet) cos I — — j cos l — sin (il$t) sm I — — I 

sm {Ugt) cos I — — J + cos (SJfltj sm \ — — 



cos L 



(8a) 
(8b) 



x 3 = r cos sin l, (8c) 

where the angular velocity in the x 3 direction is 

Vl e = M~ V [l - 12(S/M 2 )v 3 cos l] . (9) 

These equations describe a circular orbit on a plane (the orbital plane, which is inclined at 
angle i to the equatorial plane) which precesses around the spin axis of the black hole with 
angular velocity 2S/r 3 (the Lense-Thirring precession). Because the particle's motion is the 
the sum of the circular motion and the plane precession, the particle itself does not travel 
on a fixed plane, nor does it travel in a circle, nor does it cross the equatorial plane at angle 
l, but rather at an angle l' ||. 

Turn, now, to the gravitational radiation reaction force (or acceleration) that slowly 
modifies the above orbital motion. The emitted waves and their associated radiation reaction 
can be expressed in terms of the multipole moments of the "system" (particle plus hole) 
lQyilJ . In the absence of spin, S = 0, the reaction, at leading order in M/r, is due to 



the mass quadrupole moment Jy of the system; when S ^ 0, the leading order influence 
of S on the reaction is due in part to the mass quadrupole moment 7y and in part to the 
current quadrupole moment J^. The moments Jy and Jy are usually written as integrals 
over the mass and momentum densities of the source. For the black hole, however, this 
cannot be done, and even for a neutron star the standard integrals are invalid because they 
ignore the star's relativistic gravity. There is an alternative approach, however, that does 
work for this black-hole-plus-particle system: the moments are defined in terms of the weak, 
asymptotic gravitational fields far from the hole and particle. When this is done, all the 
standard formulas of the multipolar gravitational wave formalism remain valid [|1(| . 



The standard theory of the radiation reaction force (for example, Ref. ||11|| ) is generally 
formulated in the center-of-mass Cartesian coordinate system {x'j}, which differs from the 
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black- hole-centered coordinates {xj} used above. In the asymptotic, center-of-mass coor- 



dinates, the hole (or rather its asymptotically spherical gravitational field |I0 , |I2| ) moves 
along the path x' k {t) = —(p/M)xk(t) with its spin still pointing in the X3 direction, and the 
particle moves along the path x' k {t) = [1 — (fi/M)]xk{t). Correspondingly, to leading order 
in /i/M, the mass and current quadrupole moments are 

7y = \jJ,XiXj] TF , (10a) 
HXi€ jkm x k x m - Ufi/M)xiS8 j3 . (10b) 



Here, STF means "symmetrize and remove the trace", and Xi = Xi(t) is the trajectory of 
the particle in the hole-centered coordinates, as given by Eqs. (§). 

Eqs. ([10]) agree with the jj, <§C M limit of the moments given in Kidder, Will, and Wiseman 
(Ref. f7j, Eqs. (14a,c)). The first term of Jjj is the standard contribution due to the motion 
of the particle; the second term arises from the motion of the spin of the hole relative to 
the asymptotic, center-of-mass inertial frame: when the current dipole moment, due to the 
hole's spin angular momentum Jj = S5i$, is displaced from the system's center of mass by 
6x' { = —(fi/M)xi due to the orbital motion, that displacement produces a current quadrupole 
moment = [-Sx^Jj] = [second term in expression ( |10b|) 1, as one can deduce from the 



asymptotic metric components by which the moments are defined: Eq. (11.1b) of Ref. |L0 

Consider, for the moment, the radiation reaction force acting individually on the hole or 
on the particle. Let the object of interest (hole or particle) have a location x'j and velocity x'j 
in the asymptotic, center-of-mass coordinates. Then, for the moment ignoring the object's 
spin, its radiation reaction acceleration (force divided by mass) is given by 



a 



(react) _ _2j (5) / ,16 7(6)™/ , 32, I •/ , 32 t( 5 ) ^' ™' / 1 1 ) 

5 1 jk ' 45 c JPQ u pk ^q^k ' 45 c 3PQ u pk x fe x g ' 45 c pq[j u k]p • 



J 



This can be derived from Eqs. (11) and (12) of Blanchet and Damour ||11| , keeping only 
the mass quadrupole and current quadrupole moment terms. Here, the brackets denote 
antisymmetrization: Cyk] = (Qfe — Cfej)/2- The number in parentheses to the upper right 
of a multipole moment indicates taking that number of time derivatives. 
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The coupling of the object's spin angular momentum Jj to the radiation reaction field 
due to Jij produces an additional radiation reaction acceleration, 

Q>j \spin = TE^ij (^) 



where m is the mass of the object. If the spinning object were nearly Newtonian, Eq. (|P2|) 
could be derived by adding up the velocity dependent radiation reaction force [massxEq. 
flnpi on each bit of mass inside the source and then dividing by the total mass m of the 
object. For the black hole, such a procedure is invalid; however, the result flT2|) must still 
be true: The analysis of Thorne and Hartle |12|], specifically their Eq. (1.9b), shows that 



the force on any isolated spinning object in a "gravitomagnetic field" (the type of radiation 
reaction field that is responsible for this force is the same as for a spinning mass with 
only weak self-gravity. 

Since the constants of motion, E, L z , and Q, are defined and expressed in terms of black- 
hole-centered coordinates rather than center-of-mass coordinates, their evolution must be 
computed using the radiation reaction acceleration of the particle relative to the black hole, 
i.e. the difference of the accelerations of the particle and the hole relative to the center-of- 



mass: 



(react) i (react) i {~\o \ 

®j = ®j \particle O'j \black hole \LO<l) 

— ~x,Ijk X k + 4§ e jpqJpk x q x k + J^ e jpqJpk X k x q + 45 e pq[j 'Jk]p X Q Xk 15 i^/M) J^j . (13b) 

Here, the moments iy and are given by Eqs. ([[(]). There is no contribution to (|13b| ) (up to 
leading order in n) from the acceleration on the hole except through the spin interaction ( |12| ) . 
which gives the fifth term in Eq. ( |13b| ). For the in the contribution to the acceleration 
on the particle [the second, third, and fourth terms of Eq. ( |13b| )1, only the second term of 
Eq. (|10b| ) needs to be kept: the contribution from the first term is not at leading order in 
M/r for either terms not involving S or terms involving S. Similarly, the J%j in the fifth 
term of Eq. ( |13b| ) requires only the contribution from the first term of Eq. fllObp . 

By differentiating Eqs. (§) with respect to time and using the aj of Eqs. ([13]) as the 



radiation reaction contribution to Xj, we obtain the following evolution of the "constants" 
of motion: 



E = fiijdj, (14a) 
L z = ne 3jk Xja k , (14b) 
Q + L 2 Z = 2fi 2 ei j kX j x k e i i m xia m - 4 ^ Se^ k Xja k ^ 

[X m X m j 

While Eqs. (|l0|)-(ll~4]) are valid for noncircular orbits as well as circular orbits, we will 
treat only circular orbits here. Suppose that the orbit initially is circular, so the trajectory of 
the particle x k (t) in the hole-centered coordinates is given by Eqs. (§). Then, the evolution 
of its constants of motion, E, L z , and Q + L 2 , can be computed as follows. Insert the 
trajectory of the particle (|8|) into Eqs. (|T0|) to get the and moments. Then, insert 
these moments and the trajectory of the particle into Eq. (|13b| ) to get the relative radiation 



reaction acceleration. Insert this acceleration into Eqs. ( |T4"D for the time derivatives of E, 
L z , and Q + L 2 ., and average the resulting expressions over an orbit. [The averaging, denoted 
by ( ), can be taken to involve times — ir/Qg < t < +7r/Qg, and because of this restriction to 
small t, for terms in the trajectory (|3|) that have argument 2St/r 3 , the cosine can be replaced 
by 1 and the sine can be replace by its argument, thereby simplifying the calculation.] The 
result is |13 



• \ 32 u 2 in / 433 Sv 3 \ 

E) = i^v 10 1 -cost , 15a 

/ 5 M 2 \ 12 M 2 J K ' 

•\ 32 u 2 7 ( 61 - 687cos 2 ^t; 3 \ ^ „ , 

L ") = -T¥"r sl + 24 Te)' (15b) 

Q + L>) = -^ A" (l-^^ cost). (15c) 

Eqs. ( |15a| ) and ( |15b| ) agree (after trivial conversions of notation) with previous results 
(Ref. 0, Eqs. (3.13), (3.18) and Ref. [0, Eqs. (4.10), (4.11)), computed by alternative 
methods that cannot compute (Q + L 2 ). 

These time derivatives of the constants of motion have two major implications for the 
orbital evolution: First, they imply that the orbit, which was assumed initially circular, 
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remains circular; this can be seen from the fact that they preserve the circular-orbit rela- 
tionship (|7|). Second, when combined with Eqs. (U), they imply the following evolution of 
the orbital plane's inclination angle: 

ldi\ _ 244 fj,S 



v sint. (16) 



\dt/ 15 M 4 

This equation implies that radiation reaction drives the orbital plane toward antialignment 
with the spin, as might be intuitively expected since that orientation minimizes the energy 
of spin-orbit coupling |16| . 

By combining Eq. ([[(]) with the leading order change in r, r = — ^(fi/M)v 6 , the following 
differential equation relating r and i is obtained: 

dt 61 S fM\ 3/2 . 

sint. (17) 



-dlnr 48 M 2 V r J 



The time-averaging of the left hand side of Eq. ( |T7j ) is not explicitly written down. Integrat- 
ing Eq. ( |T7| ) yields 

tan(6/2) = tan(6 /2)[l + §(5/M 2 )(M/r) 3 / 2 ], (18) 

where Lq is the value of the inclination angle when the orbit has large radius. 

Eq. ( |i~8|) shows, for example, that for an 5 ^ M 2 black hole, if the particle orbits with 
a small inclination angle (0 < i tt), then i fractionally changes by 0.06 from its initial 
value at large radius to its value at r m 6M. The regime r < 6M is of special interest; it is 
there that waves from the final stages of inspiral can give high-accuracy maps of the hole's 



spacetime geometry [[17]]. However, when r w 6M, our leading order analysis breaks down, 
so that to be able to map the hole's geometry, the analysis must be carried out to higher 
order in M/r and S. 

The above analysis illustrates the power of the radiation reaction force method to reveal 
the detailed evolution of a system under radiation reaction. The case solved was sufficiently 
simple to give an easily presented solution. A future paper [|J will give more details of the 
above calculations, along with the generalization to eccentric orbits. Hopefully, future work 
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with radiation reaction forces will: (1) Generalize the analysis to an arbitrary mass ratio 
H/M and the case of both masses having spin. (2) Extend the analysis to higher order 
in M/r and in S. (3) Achieve a similar calculation of the orbital evolution in the fully 
relativistic Kerr metric. 

The author is grateful to Alan Wiseman and Kip Thorne for their advice. This work was 
supported by NSF grants AST-9114925 and AST-9417371, and by NASA grants NAGW- 
2897 and NAGW-4268. 
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